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ABSmCT 

The base condensate model of hellisa Is used to exasdne the structure of 
the electron bubble In hellua. The solution obtained makes use of the fact 
that the parameter (am/lM)^/^ la small and m/M is negligible, where m Is the 
electron mass^ M Is the boson mass» I Is the electron^^boson scattering length, 
and a Is the heallxig length* It Is shown that, to leadlx^ order, the radius 
of the bubble Is b » (rM^a^/mPo,}^^^, when Is the helium density. The 
effects of (quantum) surface tension and of polarization are discussed, and 
are shown to be small. Consideration Is given to the effective mass and radius 
of the bubble, and the elllptlclty Induced In It by slow motion Is given. The 
normal modes of pulsation of the bubble are fotmd and the mobility of the Ion 
la computed. The theory Is compared with experiment. 

1. INTRODUCTION 

It has become Increaslx^y apparent over the past decade that the dell-* 
berately Introduced Impurity can be a fruitful experimental probe of the struc-* 
ture axad dynamics of helium II, the super fluid phase of helium. Of particular 
Interest Is the tuigatlve Ion which consists of an election that, through Its 
zero point motion, ca|ves out a soft bid>ble of about 16A in radius In the 
surrounding helium (lA « lOr^cm.). The Induced hydrodynamic mass of such a 
large structure Is greatly In excess of Its physical mass, and It therefore 
responds to applied forces as would a much more massive Ion. The experimental 
situation has been reviewed by Donnelly^, and more recently by Fetter^. 

The negative Ion provides an Interesting and, as we shall see, a sensitive 
testing ground for theories of helium II. We examine In this paper one parti- 
cularly simple model of helium near absolute zero, the bose condensate.. The 
approach Is expounded by, for example, Gross^ and by Fetter and Walecka^. The 
theory Is so single to apply that most of the properties of the electron bubble 
can be calculated In an elementary way. We will present our arguments in a 
hydrodynamic framework originally proposed by Madeluag^. Since this may be 
unfamiliar to the reader. It Is developed in 1 2 for the siiiq>le single-particle 
SchrUdlnger equation. It Is generalized in 13 to the bose condensate. The two 
theories are brought together in §4, where the theory of the electron bubble 
Is developed. TbB final section (S5) confronts the theory with experiment. 

2. MADELUNGVS THANSFOHMATION 

It appears to have been Hadelung^ who first realized that SchrUdlnger's 
equation cou^d be cast Into a fluid mechanical mold, by ezpressing the wave- 
function, t}i(x,t), in terms of Its aoqilltude, f(x,t), and phase, 4(^,t). Consider 
a particle of mass m in a field of fixed potential, w^ , and therefore obeying 

ii!3V3t * -(if^/2m) wt|>. (1) 

By writing 


^ ^ i exp (Im^/fO 


( 2 ) 
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where £ and ^ are real* we can divide (X) into 


2 -|| + fV2<ji + 2^f - 0, 


■fl + 1 + a “ ^ * 0* 

When we introduce the (probabilistic) mass density » 


<3) 

(4) 


p, and current j^. 


by writing 


p . m|i|.12 - nf^, 3 - (B/2±) • ptf, (5) 

and define a velocity* u* by their ratio 

u • j/p • <6) 

we recognise that (3) and (4) are the continuity and momentum equations 
governing the potential flow (6): 


H + ^ • (p5) - G, 

If + i ^ + i + u - 0, 

where 

U - -(«^/ 2 m^) 72 p^/ 2 ^pl /2 


(7) 

( 8 ) 

(9) 


There are three main differences from classical potential flow. First, 
the total quantity of 'fluid' is not only conserved by (7): it is fixed by 

/ j tpdx * 1, or /Pdx » m. (10) 


Second, even if through the presence of walls (w»«) or otherwise the fluid 
is confined to a certain multiply-connected domain,^?, must remain single- 
. valued. It follows from (2) that, round any contour T in ^ not reducible to a 
point by a continuous deformation, can change by a snsltiple of h/m only. 

The circulation round r cannot freely take any value: it is quantized by 

the Bohr- Sommer f eld condition 

(n-0. ±1, ±2.....). (11) 

i m 

Third, a completely new term has appeared in the momentum equation (8) , 
namely y. 

The term y is often called '.the quantum pressure' • This is a misnomer 
for at least three reasons. First, its dimensionality is incorrect, and it 
would be better regarded as a chemical potential per unit mass. Second, since 


where 


6-8 


au/3x^ - 




32 

3x±3x^ 


(In 


( 12 ) 

(13) 


a rival, and properly din^nsioned, contender for the title of quantum pressure 
exists as part of the unusual and complicated stress tensor (13) • Third, the 
word 'pressure' suggests a phenomenon that depends only on 'the local thermo- 
dynamic state' (here fixed by p), and the presence of derivatives in (9), or 
(13), shorn that all neighboring points are involved in its definition. Despite 
these objections, we follow the common usage* 


As may be seen by setting if » 0 in (9) , the quantum pressure is the 
essential ingredient that distinguishes our subject from the classical 
theory. The fluid dynamicist can gather experience of its effects 
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by translating some o£ the elementary situations of quantum theory into their 
corresponding fluid iimtehanieal statements. 

* Hydrostatics* arises from the quantum mechanical bound states by writing 

^ • -Et, C14) 

i^ere £» the energy of the state, is a constant. By (8) the quantum pressure 
balances E*w everywhere. It is best to avoid the usual fluid mechanical prac- 
tice of absorbing £ into w, since some energy levels may be inaccessible. For 
example, ^en (8) is written as 

-(B^/2b) - IE • (15> 

and it is supposed that w Increases indefinitely with distance, r, from some 
origin, 0, one family of solutions to (15) is found that increase with 
r, so that the normalization integrals (10) do not converge* The condition 
that only the normalizable solutions of the other family are used transforms 
(15) into an eigenvalue problem that confines £ to discrete levels. Of course 
a continuum of eig^svalues exists when w is bounded above* 


A well-known application of (15) , that is particularly relevant to the 
bubble, is the pot^tial well for ^diich 


wx wjx 


wx* la r<b CBagioa I'); 

• la 7>b CSagloa II* ) i 

(>vx) aca onutaats. iftrltlag 


Xx2 - <2 b/«2) (B - wx). ■ <2b/«2)<wj - 


E). 


(16) 

(17) 


w« see that, for «x ^ ^ ^ obeyed by 

f - fj 5 Aji(lxr) Tl(9^X), la r < b; (18) 

f - % £ Akji(X3ir)[jt(XxJ>)/^tt<^irb)] Y^(Q,x), in r > b; (19) 


where j£(z) is the spherical Bessel function of the first kind, k£(z) is 
the modified spherical Bessel function of the second kind, and Yjt(@,x} is a 
surface harmonic of integral degree, i, in spherical coordinates (r, 9, x)* 
The exclusion of the other spherical Bessel functions ensures that (10) can 
be mat for SOTie choice of the constant A. Contimiity of f has been realized, 
and vf is continous provided 

^lMx"(^lh)/jt(Xxb) - ' (20) 

This dispersion relationship detemaines a discrete spectrum of admissible E. 
It may be seen that, when Aw = w^ - wj is large ccwapared with h^/mb^, 
eigensolutions exist for idiich E - wx « wjx E. For these, (19) takes the 
approximate form 

% w A expt-Xix (r - b)] Yit(e,x) , in r > b. (21) 

The fluid is confined in region IX to a boundary layer of thickness l/X^x , or 
- « (2mdw) (22) 
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This phenoiaenoxi is often called ^healingS the layer itself a ^healing layer’, 
and Sq^' the ’healing li^th' * To the fluid mechanicist, the abrupt increase dv 
in V at r * b can only be hydrostatically balanced by an equally abrupt increase 
in quantum pressure* It is the non-local character of quantum pressure that 
causes the fluid to pass through the barrier at r b and, when lower poten- 
tials are available externally, permits it to seep out of region I* 

In view of later developments, it is worth elaboratiz^ the situation 
just described. First we allow the well to have any shape, defining n to be 
a coordinate that measures distaime from the discontinuity normally outwards 
from I to XX* Second, we allow wx and w]x to vary, and suppose that the large 
transition f rom wj up to 'wxc does Hot occur abruptly as in (16), but contin- 
uously over some distance, a, comparable with am* Within this distance, there 
is evidently no unique way of defining ’the’ surface, S, of the well* We 
note, however, that in the case of discontinuous w just considered, jjt(Xxb) 
is small, by (20). This suggests we should locate S on the surface of zero fx* 


To elucidate the healing layer structure, we introduce a stretched 
coordinate, and cast (15) into dimensionless form by writing 

C«n/a, f-ax(5), (23) 

q • a/am, w * wx + g(5)^» (24) 

where we suppose g is exponentially small at the inner (C » -») edge of the 
boundary layer, and is unity at the outer (5 ^ +*) edge* It should be realized 
that X, wj^ Aw and g will generally depend on position on S* To the first two 
orders, however, this dependence only occurs parametrically in the solution, 
and will therefore be suppressed. Writis^ 

x ^ x^CO + ax^CC) + ..., (25) 

substituting into (15), and equating like powers of a, we obtain 

- 0 , ( 26 ) : 

+ c'Six/iU < 27 ) 

where Cxand C 2 are the principal radii of curvature of S at the point concerned* 
Since f xc is identically zero, the solutioxis to (26) and (27) must obey 

Xq 0, (28) 

Successful matching to the interior, solution, fx, requires 

*0 - I + C 2 ^^) (3fi/3n)g5^, for 5-^ — . (29) 

Explicit solutions can generally be obtained only by numerical means. They 
obey Integral conditions which we will later f±iid useful ( 14) : 


1 .afT. 2 


c H «. 


(30) 



C2^ 


■' 


dC 




( 31 ) 


where the bar through the integral sign signifies that the convergent part of 
the integral is taken* Despite appearances, this integral is negative* 
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The Integral relationships (30) and (31) may be interpreted in the light 
of the Grant tensor (13). The terms on the left give the main parts of the 
leap of <7^ across the boundary layer; the terms on the right give the 
corresponding integrated effects of the external force f^w balancing them. 

The dominant term, given on the left of (30), arises from the 'pressure* of 
the particle trying to escape the well. The n^t largest term, shown on the“~" 
left of (31), clearly has the form of a quantum surface stress, with 

\ ^ • < 32 ) 

as coeffic^t of surface tension, it is easy to make an estimate of T . If 
we take g to be a unit step function, we find from (26) that ^ 

\ <3fl/3n)|. (33) 

Passing from 'hydrostatics' to 'hydrodynamics' by abandoning (14), we 
see from (8) that another type of healing phenonmenon will occur when u is 
large. A particularly significant case occurs at a vortex line where, by (11), 
u is of order nh/mm, at small distances, m, from the vortex axis. It follows 
from (8) and (9) that p is of order S^n fo,. Jf Xhe fact that P is zero on 
the axis itself means that a closed vortex ring, or a vortex line terminating 
on boundaries (w « «),wiH transform an otherwise simply.<«oimected container 
into a multiply-connected domain,i^, so Justifying a posteriori the application 
of (11). Unlike the healing at a wall consid^ed earlier, the depression 
of p at the vortex axis occurs over distances coflg)arable with the scale of the 
container. The corresponding vortices in the condensate discussed in §3 
hive cores confined to much smaller distaxices from their axes. 


Before concluding this section we make one remark, obvious perhaps, but 
relevant to §A. ^en the particle is trapped in a potential well with moving 
walls [w * w(x,t)], ^ is necessarily non-zero and p Is time-dependent. Never- 
theless, provided the time-scales over which w changes are large compared with 
the reciprocal of the quantum frequency h/mb^, we can regard the fluid as 
being in a quasi-hydrostatic state, ignore the time derivative in (1), and 
treat t in w parametrically. In quantum language, the Bom-Oppenheimer approx- 
imation is said to apply. 

3. THE CONDENSATE MODEL 

Ve now consider an assembly of N identical particles (bosons) of mass H 
in a potential field W(^. If the particles did not interact, the wavefunction 
for the system could be written down as a aysB&etrized product of the N one- 
particle wavefunctions, f^,t), obeying (11) with W and M replacing w and m. 

It would be probably more corvenittt, however, to replace the normalization 
condition (10) by 

/|f|2^«N, or /pdx«P«v, (34) 

where v is the volume of the system and p<» « MN/v. The resulting theory is 
well-understood, and contains features that fruitfully represent helium near 
absolute zero^. It may be se^ from (18), however, that the ground state for, 
say, the potential well (16) would be one in which all the particles would 
be at the origin, with high probability. To eliminate this unphysical behavior, 
the imperfect bose condensate has been devised. A short-range repulsive 
potential V(? - is introduced in an ad hoc way, and to W the potential 

/V(? - t') I ?^') 1 (35) 
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Is added which Increases as the density of nearby bosons increases. The 
simplest case arises when V is taken to be 

V(x - x') - «(x - x'). (36) 

"Equation (X) is then replaced by 

ai3f/3t - -(-Ii2/2M)v2¥ + C37) 

A fuller and oore satisfactory derivation of (37) nay be found elaewhere^'^. 
It is of some interest that non-linear Schrodinger equations of the form (37) 
have been the subject of close scrutiny in recent years in non-quantal 
contexts, particularly in theories of weak non-linear waves and. stability^. 
The Hadelung transformation 

Y « Fexp(iM^/h), (38) 

follows the course of §2 with minor changes. Host significant is the 
addition of a * gas pressure^ , 

p - (Vo/2M2)p2, (39) 

which (multiplied by 6 ^j) should be included in the stress tensor (13) * 

Thus (8) is replaced by 

H + i ^ + i + ■ 0- 

The presence of the repulsion, and its associated gas pressure 
restores a number of physical effects absent in §2. The tendency towards 
condensation is eliminated for all sufficiently large systems. To see 
this, return to the hydrostatic theory of 12 and the potential well (16). 

The spherically symmetric (i«0) grotmd state now obeys 

- if ■ 

If Wq and AVf are both large compared with (41) gives everywhere 

except near the surface, S, of the well 

p » p« « MFx^ « MN/v, (42) 

the corresponding one-particle energy being given by 

E - Wx - P^Vq/M. (43) 

The fluid is spread out uniformly in the well. 

Near S the derivatives of F become large, and the constant solution 
(42) breaks down. We may follow the argument of §2. Introducing a new 
healing length 

a - ^44) 

writing 

C - n/a, F - x(5), W - + (p,Vo/M)G(5), (4S) 

where 6(0 is exponentially small for Z-*- - *>, eicpanding Z as 

X(0 - X^(0 + aXi(0 + ... . (46) 

substituting into (15) > and equating like powers of a, we obtain 
d2Xo/d52 - [6(0 - 1 + (3 Co/Fi)2 3 X^ - 0, 


( 47 ) 
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d2xi/dc2 - [G(0 - 1 + 3 (Xo/Fj) 2] - - (Ci"! + C2"^)dXg/d5. (48) 

Hatching at tha edges of the boundary layer requires 

?X* H - •. (50) 

Again, explicit solutions generally require numerical integrations, though 
useftil integral relations may be established, for Instance 

i - -- f « • <5i) 

-2<i +i> il « • iZ 2X^^d5. (52) 

1 2 dC 

Once more, the dominant contribution to the leap in stress across the healing 
layer arises from the interior solution, although now it is the gas pressure 
and not the quantum pressure that is mainly responsible* Again, the next 
largest term can be interpreted as a quantum surface tension, with positive 
coefficient 

Xm - 2a(£^) r (^)2 dc . (53) 

As before, % may be estimated^ from a simple model of X • If we suppose 
that G is a step function of infinite height (AW « •) , (?6) may be solved 
as Xq » ?x (rV^t) 9 and (53) gives 

i /2 -K^p (54) 

Passizig again from * hydrostatics' to 'hydrodynamics' , we note that 
the gas pressure can supply the restoring force necessary for coaq>ressional 
waves* Perturbing about the static solution (41), we readily find that 
long wave-length sound propagates at the velbcity 

c • /(dp/dp) * /(2p/p) » -h/Ma/2* (55) 

At wavelengths of order a and smaller, the quantum pressure increases the 
phase speed, decreases the group velocity, and introduces weak dispersion* 

Vortex lines may be* studied as in §2* Unlike their classical counter- 
parts, the cores of these vortices do not have sharply defined surfaces 
separating regions of zero axui non-zero vorticity* All the vortlcity 
they contain is concentrated as d-functlons on their axes* Such a vortex, 
if classical, would have infinite self-energy. Here, however, the density 
decreases over the characteristic distance a as the axis is approached, 
so ensuring a finite tension. the depletion of fluid in the core makes the 
vortex resemble the classical hollow core model* the permanence of vort» 
rings implied by the Kelvlxb-Helmholtz theorem makes them excellent candidates 
for quasi-particle models, so reviving in a novel context the ideas under- 
lying the vortex atoms proposed by Kelvin in the nineteenth century. 

4. STRUCTURE OF THE NEGATIVE ION 

It is possible^ to account with relative ease for many features 

of the negative ion by combining the methods of l§2 and 3 above. We use 
the theory of §1 to represent the electron, regarding w as the potential 
created by the surrounding condensate; we apply the formalism of §2 to the 
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exterior of the bubble, taking for W the, potential of the electron. Hore 
explicitly, we introduce the energy, 

//d( 5 - r> |<i(5E) I (56) 

representing the repulsion of an electron at x and a boson at Taking 
again the sinplest case of a 6<»function interaction 


u(t - r) - ti^6(* - r). 

(57) 

we then have^ ^ ^ 


w(x) - /U(S - xOlT^') |^dx' - 

(58) 

wet') - mt - 1') I'l'et) pdx - Uoiii<(r)|2. 

(59) 

By (10) and (58) , Aw « U©p«/M so that by (22), (23) and (44), 
Equations (1) and (37) become coupled: 

• vSJg/mp. 

Ur8iK/9t - -(«2/2m)72t|. + Uo1t|^. 

(60) 

OiiVit - -(«2/2M)v2y + (Vo|f|2 + Uo1'('|*)T, 

(61) 


As in §2, we define the electronic surface, S, of the bubble by the zero of ^ 


The key to a simple * hydrostatic^ solution of (60) and (61) lies in the 
fact, which we can verify a posteriori, that the radius, b, of the bubble 
is large compared with a and so that the boimdary layer methods of §§2 
and 3 can be used with minor emendations. We must not forget however that, 
since the roles of Interior and exterior of the bubble have been exchanged 
for the condensate, the sign of C in §3 must be reversed. The mainstream 
value of F, denoted in 13 by Fj is now written Fg. 


To leading order, we set g(5) * X^^/Fs^ and 0(C) • 

(26) and (47), as (38) and (59) require* The integral relations (30) and 
(51) may then be combined to give 


2 


fii) 2 




- 0 


(62) 


where we have appealed to (28) and (50). To the next order, the forms of 
g and 6 require reconsideration^^. In place of (27) and (48) we have 

d2xi/dC2 . (q/Fs)2(Xo^ + - -(Cj^"^ + C2”^)dXo/d5, (63) 

d2Xi/d€2 - OXo^/Fg^ - DXj, - (a2Do/Fs^Vo)(xo2Xi + 2X0X0X2) - 

- -(Ci-2 + C2"^)dXo/d5. (64) 


The integral consequences (31) and (52) are modified accordingly, and the 
result (62) is altered to 


«2 


(|£) ^ - -2o 

'^an^S 2W„ 



which now includes the effects of interfacial tension. 


The jusp conditions (62) or (65) across the boundary layer suffice to 
match the mainstream electron solution to the mainstream condensate solution. 
Applied to the electron bubble, we have by (18) 

f ■ A(sin ^x^)/^X^, 


( 66 ) 



43 


ffaowias that Xx^ * ^ solution is normalised to the first two 

orders^ in a/b, if * ir/2b3. ^ (10) and (U) we obtain 

^ ^ «• <« 7 ) 

It is not at once clear idiether b will be decreased by the positive surface 
tension (S3) of the condensate or increased the negative swface tension 
(32) of the electron. If we use estimates (33) and (54) however, we see that 
1^1/% is <>^der <C^ and, since erperla^nts indicate (§5) that g < 1, it 
appears that the bubble radius shouixi be larger than ('nH^a^/mPw)^'^. 

Direct numerical integrations^^ of Xq and and evaluation of the integral 
seen in (67), suggest that the dtf ferries is of the order of a. 

The effects of polarization induced by the electron in the surrounding 
helium can be included by addingi^ to (56) the term 

|2|5f , ?-j”^dxdx', (68) 

where a is the polarizability of the helium and e is the electronic charge. 

This has the effect of contracting the bubble by order o4!e2a2/4irR^b3. A 
detailed theory^^ shows that the reduction is of order a/3 in the practically 
interesting cases* 

Further complications arise ^en the dynamics of the bubble are considered, 
although the time-scales of interest are usually large enough compared with 
the electronic frequencies to justify the neglect of in (60); see 12. 

To evaluate the effective hydrodynamic mass, we consider the bubble in steady 
motion 0, at small Mach makers ^0/c. The electronic radius, b, of the 
bubble is increased by about 5^Z because of the pressure forces associated 
with the flow of condensate over its surface* It is also made slightly oblate, 
wltir an. elliptlcity close to S^/2. Ignoring this effect, it is found that the 
dipolar back-flow created by the ion coincides with that of a hard sphere 
whose radius, be» is less than the electronic radius, b, by one to two 
healing lex^ths: 

b* - b - <aM/p-)^Xo2(5)d5. (69) 

It is this radius, rather than b, that determines the induced mass of the ion. 

Further details of the calculations mtlined above may be found in the 
paper by Roberts and Grant already cited^^. We conclude this section by 
breaking new ground* We consider the oscillations of the bubble, their impli- 
cations for phonon-ion collisions, and tim mobility of ne^tive ions at low 
temperatures* We again adopt the boundary layer methods described above but, 
of course, retain t^e term 3Y/3t in (61), so introducing a velocity potential, 
in the condensate. We retain only the dominant part of the boundary layer 
structure, excluding both surface tension and polarization effects* We write, 

f » f ^ 4- oaf", F • F^ 4- oaF", ♦ • oO^, (70) 

where the suffix o stands for the steady solution obtained earlier, and the 
terms in a represent time-dependent perturbations ,idiere 0 < a « 1* 

It is readily seen from (7) and (40) that and F^ both obey the acoustic 
wave equation 

326'/3t2 . 


(71) 
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13 

We could » by foUovlng Celli, Cobeu and Zuckeraian , examine solutions In tbe 
form of outgoing mves. the elgenfrequencles would be cosqtlen* because of 
the reduction In oscillation amplitude at a point fixed In space as the 
energy of surface motion Is radiated to Infinity, the corresponding elgen^ 
function must tend to Infinity with r slnce» the more distant a wave ls» 
the earlier It must have left the surface, and the greater the amplitude 
of surface oscillation most then have been* We will not consider solutions 
of this type below* We will confine our attention to the scattering problem 
In which an Incoming plane wave travelling In the z -direction 

- eai»Ii(kz- a)t)l, (72) 

where Ic and m * ck are real. Is scattered by the bubble Into a set of outgoing 
waves. We first aim to calculate the scattering asq>lltude, hji(k), of the 
1th partial wave*^ 




1 

1 + IQn* 


Qt 


kby»(kb) + K(»vii‘‘(kb) 
kbjji(kb) + Kjjji'(Kb)’ 


(73) 


where 7j^(z) is the spherical Bessel function of the second kind and is 
the spring constant of the bubble for this mode. We then use bf, to coi^ute 
the differential cross-section^^ of the bubble 


o(k,0) - k"^| Jq( 21 + l)h)i(k)P^(cose) P, 
and hence the momentusk-^transfer cross^sectlon^^ 
of^(k) » //<^(k,0)(l - co80)sln6d0dx» 

From this we finally evaluate the mobility, Ue» fromf 

where 


14 


n(k) - [exp(Kck/KT) - 1 

is the density of phonons (in ^^space) at temperature T. Here K Is 
Boltzmann's constant* . 


(74) 

(75) 

(76) 

(77) 


To determine the spring constants, we have to match solutions of 
(71) across the boundary layer on S', the deformed electronic surface, to 
the quasl^-statlc solutions (18) of the electron mainstream. We first consider 
the case t > 1* The fact that 

* - AoJodir) + txaA'iiazT) Yji(0,x). (78) 

Implies that has the equation 

r - b(0,x.t) 5 b^ + aab"(t)Y^(0,x), (79) 

where ' (using X jb^ « v) 

^>"/^o - itW AVAo* (80) 

We will continue to refer the boundary layer structure 'to the unperturbed 
position, Sq, of the electronic surface, and not to S^* We Introduce x', X' 
and n the boundary layer forms of f, F' and and expand these in 
ascending powers of a 

X * a Xq + Xj^ + ••*, X “a Xq 4* Xj^ H* ••*, ^ 

( 81 ) 
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where the coefficients shown depend on t wad I (r - ho) /a and paraTOtricslXy 
on @ and x* We substitute these into the boundary layer forms of (60) and (61) 
which, in the Mhdelung framework, give to the first two orders 

- -(2a/b) 9x'/ 3g, (82) 

32xV35^- OXoVFs^ - 1)X' - (a2VFg2^|(3^2x' + 2Xo*oxO^^- 

- -(2a/b)3X7 35 + (M/Fs2Vo)Xo3n'/3t, (83) 


Xga^Ti'/ag^ + 2(dXo/d€)3n735 - -(2a/b)Xo3n73g - 2a^3X73t. (84) 

It seems clear from (84) that ^^3no73g la independent of S and, since 
there is no net eond^isate flow through the boundary layer at any point, that 
constant must be zeroe Thus Po takes throui^ut the boundary layer the main- 
stream value $s of evaluated on Sp. right-hand side of (83) does not 
contribute to leading order in a, and (82) and (83) may be solved to give 

xb" - Cd3^/dC, Xo" * ?dXo/dC, (85) 

where t is Independent of These forms represent a net displacement of the 
equilibrium boundary layer from Sp to S^, without change of form; we conclude 
that t » -b"tjf 

In proceedii^ to the neact order we note that, since the velocity of 
sound (55) is of order 1/a, the time derivatives in (83) and (84) now 
contribute. In fact, excludix^ again a net flux of condensate through the 
boundary layer, (84) shows that 3ni^/3€ takes the value -a^dc/8t throughout, 
and in particular on the outer edge (€ « «*) of the boundary layer. It 
follows that 

a3b73t - (3478r)s, (86) 

an, equation with an obvious interpretation. The equations (82) and (83) 
again admit an integral, namely 


On taking the limit C and using (18) to evaluate the contributions from 

the lower limits, we find 

34g'/3t - c^Kjiab'/bot (88) 

idiere 


f 5/2, if 

U -4 (89) 

( 1-i + if 4 > 1. 

The numerical values of Yii for the first 20 values of i are given in table I. That 
of Xp was obtained from an analysis too similar in spirit to the one just 
described to be repeated here. It may be noted that is zero, representing 
the fact that the bubble is neutrally stable to a uniform displacement. 

Equations (86) and (88) are applied on Sp, and provide the boundary conditions 
to which solutions of (71) must be subjected. 


We developed a program for an Hewlett-Packard 9820 A desk cong>uter to 
evaluate o^(k) and Ue from an arbitrary set of the spring constants* The 
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results were tested with the values (K^ * 0.23474; % * 0; K 2 * 0.45045; 

« 0, u>2) used by Baym» Barrera and Pethlck^^, and good agreen^nt was 
obtained. The prograome was then used to generate the results shown in 
Figures 1 and 2. The effect of truxu^ating s eries ( 74) at i «» 2 and £ « — 

are shown in both cases. The prominent peak in <Jx(k) seen in Figure 1 is 
due to a d-wave resonance (i » 2). A new minor peak is added every time i 
is incremented by X. The curve appears to approach its geometrical value^^ 
(0x/4vbQ^ 1/2) with an oscillation of amplitude (kb^j)*2/3 and period 

In Figure 2, we see plotted in units of Lq, as a function of T measured in 
units of To where 


To 

- «c/boK, Xo 

« 3ireis2c3/2bQK. 





TABLE 1 





Spring' Constants 



1 

U 


t 

% 

0 

2.500000000 


10 

U. 56327795 

1 

0.000000000 


11 

12.59928X01 

2 

2.289868134 


12 

13.62973970 

3 

3.771253431 


13 

14.65585350 

4 

5.032253885 


14 

15.67849696 

5 

6.198547165 


15 

16.69832328 

6 

7.314641577 


16 

17.71583055 

7 

8.400646541 


17 

18.73140537 

8 

9.467085072 


18 

19.74535255 

9 

10.520037400 


19 

20.75791573 


5. FXPERB!£!n:AL COMPARISONS 

The condensate sg)del of helim II is essentially a theory having only one 
disposable parameter, namely the pseudo-potential, V , or equivalently the 
healing length, a. It is natural to seek to choose this so that theory and 
observation are in optimum accord. Clearly a choice of a made to fit one 
physical phemomenon well is likely to conflict with others, and an overall 
consistency with the experimental facts is not to be anticipated. One notes 
particularly that, since the condensate is a gas obeying the equation of state 
(39), we should not expect the theory to perform well at the vapor pressure. 

One can obtain an estimate of V ■ 4trd!i^/M from measurements of the atomic 
diameter, d, by a-particle scatteriz^ experiments. Values of d of about 2.7 A 
have been found. If is 0.145 g/cm3, the healing length would be 0.82 A, 
leading to too small a velocity of sound. One popular procedure has been to 
extract a from accurate experimental determinations of the relation between the 
velocity and energy of circular vortex rings. This had led to estimates of 
a *t 1. 28 A, giving much too small a value of c. The reliability of the approach 
can, however, be questioned. One would have expected a to decrease with 
increasing pressure, but the reverse seems to be true^^. It is now believed^7 
that the core of a superfluid vortex is the seat of excitations (normal fluid), 
and that the surface of a vortex core marks the distance from the axis at which 
the Landau critical velocity is reached, rather than a quantum healing distance; 
such a belief is consistent with the increase of a with pe». 
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Perhaps the oiost satisfactory way of estlioatixig a Is through the velocl^ 
of soui^ (55)* to give an eran^le, if we take c * 238 a/a> we obtain a f 0.47 k 
and for po» » 0*145 g/ot^ ^ that % * 1.7 10*37 g o^/a2» a va3m 
admittedly three times larger than the scattering experiments suggest* and 
moreover one which will alter as c and p« change through applied pressure* 
Neverthelesst by using c to determine a» we obtain a coefficient of surface 
tei^ioa, %, from estimate (54) of 0*37 g/s^* ’«hich is in good agreement with 
the eaqperimeatal value of 0*34 g/s^ at low temperatures* 

turning to the bid)ble* we see that* in the first app^ximation* tlm theory 
does not require a knowledge of the pseodo«potential, XTo, either for its equi- 
librium structure or for its oscillation spectrum, the radius* (trH3a3/mpo»)l/5, 
predicted by the first approximation is somewhat small* 11*8 & using the viilues 
qwted above* Since c and p» increase with ixuereasing pressure, p* this radius 
decreases with increasing pressure* although somewhat more slowly than ^eri- 
ments indicate* the bubble radius is increased when the effects of surface 
tension are added* Unfortunately, Grant and Eoberts^^ did not examine values 
of a as small as 0*47 1, so that the value of the integral appearing in (67) 
is not known* Using our^earlier estimates, however, it appears that b will be 
increased to about 13 *3 A by surface tension effects* table 2 gives u^t^ for 
a few values of t for both the 1*2 and the 1 * 19 truncations, and for values 
of b of 11*8 I, 13*3 1 and 16*0 A* At the 1*2 level of truncation, there is 
a clear tendency for PeX3 to approach a limiting value, of about 36 cm^K^Vs in 
the case of the 16 A bubble, as t increases* the explanation of this behavior 
was provided by Baym, Barrera and Pethick^^ in terms of the shape of the d<«wave 
resonance of Pig* 1* Hot surprisingly in view of the very different form of 
Of obtained at 1 • 19 truncation, the constancy of y^t^ is not as marked at 
this level* Heven^eless, the values shown for 1 * 19 In table 1 are not ridi- 
culously far from the experimental value^3 of about 32.5 cm3K3Vs in the range 
of t in which Baym» Barrera and Pethick measured the success of the work* 


When we take the th^ry of the bubble to the second approximation, a new 
disposable parameter enters, namely the psuedo-potential, Uq, or equivalently 
q » a/ajQ, a relation we can also write as q^ * mUo/M^o • Roughly speaking q, 
as the ratio of the two healing lexigths, measures the relative penetration of 
the condensate wavefunction into the cavity to the penetration of the electron 
wavefunction into the condensate* If q were zero, it would be legitimate to 
treat the condensate as an abrupt edge and only consider the electronic boundary 
layer of §2* At first sight it might appear that, since q^ is proportional to 
m/H f Io37 icr^ it would be admissible to follow Celli, Cohen and Zuckexman^3 
in taking this view, the indications are, however, that Ug/Vo is large* 
Scatterix^ experiments give im electrox»-»helium scattering len^h. A, of about 
0.60 A, iaqilying that U© * 2if>S^/m±s about 4.6 10*35 g cmS/s^. taken with 
the ea^erimental value of 5*7 10“3® g cafi/s^ for we obtain Uq/Vq A 810 
and q a 0*33* It would be intaresting to see idiether the effect of restoring 
q to the CeUi-Cohen-Zuekerman theory would have serious repercussions* the 
indications are that it would not* 


the neglect of q in the condensate theory described here would eliminate 
the condensate surface tension, and transfom the interfacial boundary 
layer into the structure considered in §2* the associated negative surface 
tension', t^, would tend to expand Che bubble, an effect confirmed by the cal- 
culations of Grant and Boberts^^. influences of Interfacial tension are, 

however, of second order in the condensate theory* In the approach of Celli, 
Cohen and Zuckerman3‘3, the interfacial t^islon, a, la a first order effect* 
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These authors regarded o as a disposable parameter that eould be legitifflately 
chosen to fit the obsenred bubble radius at the applied pressure of interest* 

It is easy to irerify that the interfaeial tension they require is positiire 
and« particularly at higher pressures, sefueral times larger than the conden-* 
sate surface tension, 1^, considered earlier. It may be wondered mhy, with 
this sign difference, the bubble radii th^ obtain are, being in perfect 
agreement with experiment, larger than those obtained from the condensate. 

The answer is to be found in (39). The condensate is a gas and, to obtain 
agreement with the observed helitmi densities, it is necessary to choose a 
large 7^, leading to pressures (39) of the order of 40 atmospteres. con- 
trast, the Celli, Cohen and Zuckerman theory treats the helium as a classical 
cooq^ressible fluid, not containing the pressure (39), and to avoid large 
bubble radii at the vapor pressure, a positive interfacial tension is needed. 

As we have stated above, we regard (39) as an artificial construct of the 
theory, not to be identified with the applied pressure, and base our comparison 
with experiment on density and velocity of sound data. 

As Schwara^^ observes, if the spring constants were regarded as disposable 
parameters, there would be no difficulty in reproducing any ion mobility data 
precisely. It appears that even the added fle^bility given to the theory by 
the ad hoc interfaeial tension, o, already permits an excellent account of the 
mobilities. Schwarsi^ has shown that, for their spring constants, the constancy 
of in the range of T of interest is not lost idien the^truncation level is 
increased as it is in ours. In comparing our theory with theirs, one must be 
perplexed by the substantial difference in the spring constants and in the shape 
of the mobility curve (labelled 19' in Fig. 2). He must wonder if , in the 
disappointing form of that curve on the present theory, and in the sensitivity 
of the mobility itself to the healing length {es evinced by the e-dependence 
of (90)1, the condensate theory has not met its most severe test to date. He 
may also speculate on the physical basis of the ad hoc interfaeial tension re^ 
quired by the other approach to survive its trial by experiment, and also 
whether the effects of roton-ion collisions at the higher temperatures have 
been underestimated. 


TABLE 2 

Ion Mobilities, 

(T in degrees K, y^T^ in units of 



X » 2 truncation 




1 * 

19 Truncation 

bo-16 

.0 1 

bo-11 

.8 A 

bo-13 

.3 A 

bo-16 

.0 A^ 

bo-U 

.8 A 

bo-13 

.3 A 

T 


T 

y^T^ 

T 


t 


T 


T 


0.34 

229. 

* 0.37 

117. 

0.31 

97.0 

0.32 

265. 

0.33 

158. 

0.35 

70.5 

0.41 

131. 

0.45 

77.2 

0.43 

51.8' 

0.37 

178. 

0.42 

84.9 

0.39 

59.0 

0.51 

86.7 

0.52 

62.3 

0.53 

41.7 

0.48 

95.3 

0.46 

71.0 

0.43 

49.0 

0.58 

70.0 

0.56 

56.9 

0.58 

39.5 

0.52 

79.7 

0.52 

59.0 

0.50 

40.0 

0.63 

63.9 

0.64 

50.2 - 

0.65 

37.6- 

0.58 

66. Z 

0.57 

51.1 

0.59 

33.1 

0.72 

56.3 

0.70 

47.5 

0.73 

36.6 

0.64 

57.3 

0.71 

39.8 

0.70 

28.1 

0.78 

53.4 

0.78 

45.2 

0.78 

36.5 

0.80 

44.7 

0.85 

33.8 

0.83 

24.6 

0.98 

49.5 

0.94 

43.9 

0.95 

37.4 

1.01 

36.0 

1.01 

29.6 

1.02 

21.0 

1.18 

49.7 

1.15 

45.0 

1.16 

39.9 

1.13 

30.1 

1.13 

26.8 

1.14 

19.4 
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